Non-Markovian Effects on the Geometric Phase 
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The geometric phases of a two-level atom interacting with non-Markovian environments are calcu- 
lated and the non-Markovian effects on the geometric phases are discussed in this paper. Three kinds 
of methods that describe the non-Markovian process, projection superoperator technique, memory 
kernel master equation and post-Markovian master equation, are used in the discussions. The re- 
sults show that when the dissipation rate is large, the non-Markovian effects change the geometric 
phase more strikingly than the small one. 
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I. INTRODUCTION 

The so-called Berryphase or geometric phase was first 
pointed out by BerryjlJ when he studied a pure quantum 
system that undergoes an adiabatic cyclic evolution. Si- 
mon gave an interpretation of this phase in the language 
of differential geometry and fibre bundles 0- Since then, 
large numbers of work have been done to generalize the 
conception of Berry phase. For example, Aharonov and 
Anandan extended the Berry phase to the case of non- 
adiabatic evolution Samuel and Bhandari [1] general- 
ized the Berry phase in cyclic evolution to the case of 
non-cyclic evolution. 

Another direction of the generalization is to study the 
geometric phase for mixed states and nonunitary evo- 
lution. This work was first done by Uhlmann[5[ via the 
mathematical context of purification. Sjoqvist et at. gave 
an alternate definition of geometric phase for nondegener- 
ate density operators undergoing unitary evolution based 
on quantum interferometry|6[. This study has been ex- 
tended to the case of degenerate density operators by 
Singh et al. 0] ■ Geometric phase in nonunitary evolution 
have been addressed in Ref . [H, Q . Recently a kinematic 
approach to the geometric phase for mixed quantal states 
in nonunitary evolution is proposed by Tong et aZ.[To|. 
The conception of off-diagonal geometric phase for pure 
states was proposed in Ref. [Ill] and has been extended 
to mixed states in Ref. [12]. The relation of geometric 
phase between entangled system and its subsystem has 
been studied in Ref. [13] . 

The intrinsic nature of geometric phase provides 



us an impelling tool for fault tolerance quantum 
computationfTij . However, the unavoided interaction be- 
tween system and environment would destroy the co- 
herence of the quantum system and hence to limit the 
implementation of quantum computation. So the study 
of environment effects on geometric phase is highly re- 
quired. Geometric phases in dephasing system and spin 
environment were studied in Ref. (l5j , which used the ex- 
act solution of the system dynamics and then it is fully 
non-Markovian. However the exactly solvable models for 
open system are only few, most problems are treated with 
certain approximations. For instance, the finite temper- 
ature effect of environment on the mixed state geometric 
phases [l6| and the effects of a squeezed vacuum reservoir 
on geometric phases [13] were studied within the Marko- 
vian approximation. It is well known that when the dissi- 
pation is large or reservoir is finite and the initial state of 
the composite system(system and environment) is entan- 
gled, the process exhibits strong non-Markovian effects. 
In this paper we shall study the non-Markovian effects on 
the geometric phase of a two-level atom interacting with 
a non-Markovian environment. Three kinds of methods 
which describe the non-Markovian process, projection su- 
peroperator technique, memory kernel master equation 
and post-Markovian master equation, are considered. We 
calculate the geometric phase for the system governed by 
these equations and discuss the results obtained with dif- 
ferent dissipation rates and different memory effects. 

Throughout this paper, the geometric phases are cal- 
culated according to the formula of TongjlOj. 
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where Xi (t) are instantaneous eigenvalues of the density The dot on ipi (t) denotes the derivative with respect to 
matrix p(t) and \ip.i(t)) are corresponding eigenstates. time t, Arg expresses the argument of a complex number. 
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of energy levels which are arranged in two energy bands 
with the same width(see Fig[T]). The levels of each band 
are equidistant. The lower energy band contains Ni lev- 
els while the upper band N 2 levels. The transition of the 
atom is in resonance with distance uj between the bands 
(we set % = 1). The total Hamiltonian for such a system 
is HU H = H s + H E + V, with 
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FIG. 1: A two-level atom coupled to an environment consist- 
ing of two energy bands with a finite number of levels. 



V = A 2J c(ni, n 2 )cr + \ni) (n 2 \ +H.C. 



(3) 



The structure of this paper is organized as follows. In 
Sec. II we discuss a two-level atom interacting with a two- 
band environment. The evolution equation is obtained 
via the projection superoperator methods. In Sec. Ill, we 
calculate the geometric phase in an open system governed 
by the memory kernel master equation and in Sec. IV gov- 
erned by the post-Markovian master equation. Conclu- 
sions and discussions are presented in Sec.V 



II. NON-MARKOVIAN EFFECTS WITH 
GENERALIZED LINDBLAD MASTER 
EQUATION 

In this section, we will calculate the geometric phase of 
a two-level atom(system) coupled to a two-band environ- 
ment in both Markovian and non-Markovian cases. The 
evolutions of the reduced density matrix for the atom 
in both cases are obtained through correlated projection 
superoperators method [Til fl9l I2Q. l2ll I22I] . different cho- 
sen projection operators would differentiate the Marko- 
vian and non-Mardovian cases. To calculate the geo- 
metric phase we solve the time evolution of the system 
in the interaction picture then back to the Schrodinger 
picture according to the free Hamiltonian of the system 
Hs = \wa z . In the following discussion, the initial state 
is chosen as, 
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Where |1) and |0) are the excited state and ground state 
of the atom respectively. This state corresponds to a 
state vector in Bloch sphere with polar angle 9 and az- 
imuthal angle cb. 



A. Model 

Consider a two-level atom coupled to a two-band envi- 
ronment E. The environment consists of a large number 



where the index n± denotes the levels of lower energy 
band and n 2 the levels of upper band, a z and a + are Pauli 
operators, A is the overall strength of the interaction, 
c(ni,n 2 ) are coupling constants which are independent 
to each other and satisfying, 



(c(ni,n 2 )) = 0, 
(c(n 1 ,n 2 )c(n' 1 ,n' 2 )) = 0, 
{c(n 1 ,n 2 )c*(n[,n' 2 )) = $ nu n'5n 2 ,n' 



(4) 



Setting Hq = Hs + He, we transform our discussion 
to the interaction picture that, 



V(t) = a+B(t) +a~B\t), 



(5) 



where 



B(t) =\J2 c(n 1 ,n 2 )e- to (™ 1 ' n2 )*|n 1 )(n 2 |, (6) 



and 



w(ni,n 2 ) = Se 



n 2 m 
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B. Markovian case 

When the standard projection is used, the evolution 
of the reduce system must be Markovian (lH. H9I I20I ]. If 
we project the composite system's state into the form 
Vp = {Ttep)®Pe — Ps® -j^-IIi, where V is the projection 
superoperator and IT = J2 n± |ni)(m|, the second-order 
Markovian master equation in the interaction picture is 
obtained [13| : 

j t Ps(t) = 72 (a-p s (t)a+ - a' , p s {t)}\ , (8) 



where 
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This equation is the same as that describes a two-level 
atom spontaneously decay in vacuum under the Marko- 
vian approximation. The time dependent of the reduce 
system in the Schrodinger picture may be obtained easily 
by a simple calculation. With the initial state Eq.([2]), the 
final state follows, 



It is obviously that the eigenvalue A_ = at t — 0. From 
the formation of Eq.([T]), we can see that the eigenvalue 
A_ and its corresponding eigenstate |— ) contribute zero 
to the geometric phase. This simplifies the calculation. 
The eigenstate correspond to A+ can be written as: 
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The two eigenvalues of pit) are 
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where 
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where 
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Obviously, for t = 0, tan If = cot|, as expected. Now 
we substitute Eas- fTU]) and (fT5|) into Eq.([T]) to obtain the 
geometric phase at time t, 
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We consider the geometric phase acquired after a 
quasi-periods T — 22. In this case, we can rewrite the 
geometric phase as 



<!>,;/. : / COS 2 -±dt, 



where cos 2 If can be obtained according to Eq. |13 
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Non-Markovian case 



If we choose a correlated projection on the total system 

Vp = Tr B (nip)®-^ni + Tr s (n 2 p)®— n 2 

iVl iV2 
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where 1I2 = l n 2}("-2|, the evolution of the reduced 
system is said to be non-Markovian \l§L [l9l . [20j . At this 
stage the master equation in the interaction picture is[22j 



j/s\t) = l2V- P [ s\t)a + - 2L{a-a+, p f(t)}, (17) 



where 71 has the similar definition to 72 and the state 
of atom is given by p$ = p)p + pg\ This equation is 
of a generalized Lindblad form[l9j and gives an excellent 
approximation of the reduced system's dynamics [22|. 
We assume that initially only the lower band is popu- 

(2) 

lated. This means p s (0) = 0. We also set 71 = 72 = 7. 
The evolution of the atom with the initial condition $2$ 
in the Schrodinger picture is 



p(t) 
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Following the same procedure mentioned in Marko- 
vian process, the geometric phase acquired after a quasi- 
periods T — — takes the same form as Eq. (fT5|) . where 
9t and rj are defined by 
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The geometric phases as a function of initial polar an- 
gle 9 in both Markovian and non-Markovian cases are 
shown in FigO In this figure, we have set transition 
frequency u> = 1 and 71 =72 =7- The geometric 
phase <&gp an d the polar angle 9 are plotted in units 
of 7r. We see from this figure that when the dissipa- 
tion is weak (7 = 0.1), the geometric phase for both pro- 
cesses are similar especially when 9 > 5, whereas when 
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the geometric phase is usually interpreted as the solid 
angle of the evolution track in Bloch sphere. We also 
note that in the figures here and below in this paper, 
the left boundaries of 8 are all set 8 ~ * 0. If 6 = 0, 
the geometric phase in all these processes are zero. The 
qualitative analysis (see the solid line in the bottom fig- 
ure) is following. In the Markovian case, because of the 
semi-positivity of ry, the two branches of the eigenvalues 
A± maintain the relation A + > A_ all the time for 6^0. 
The calculation of the geometric phase given above is 
available. When 8^0 and jt — In 2 the two eigenval- 
ues tend to be equal, i.e. the eigenstates of the density 
matrix are approximately degenerate. After this critical 
point, the relation A + > A_ holds again in the above 
calculation. However, as a matter of fact, in the case of 
8 = 0, the two branches should be crossing, and the rela- 
tion A+ > A_ should be changed after the critical point. 
Moreover, when the degeneration occurs, the geometric 
phase should be calculated as discussed at the end of 
Ref. [nj, which in our model shows that when 8 = 0, the 
matrix elements of p(t) and the eigenvalues and as well 
as the coefficients of the eigenstates are all real. Those 
facts make the geometric phase zero all the time. This 
results in the difference in the two cases 8 = and 8 — > 0. 



III. NON-MARKOVIAN EFFECTS WITH 
EXPONENTIAL MEMORY 



FIG. 2: Geometric phase in a system governed by the gener- 
alized Lindblad master equation. The geometric phase $gp 
and the polar angle 9 are plotted in units of it. We set the 
transition frequency uj = 1 and time r = — . The parame- 
ter 7 is chosen as: top: 7 = 0.1, bottom: 7 = 1. The solid 
line indicates the Markovian process while the dash-dot line 
denotes the non-Markovian one. 



the dissipation is large(7 = 1), the non-Markovian effect 
changes the geometric phase drastically, especially when 
6^0. Because the time dependent for the off-diagonal 
elements p±2 and P21 are the same, we shall study the 
difference of diagonal elements p\\ in both cases. The 
decay functions are Am = e~ 7t and An = |(1 + e -27 *) 
for Markovian and non-Markovian equation respectively. 
When jt is small enough, we can make a Taylor expan- 
sion and ignore all items higher than the second-order. 
Then Am = Ajv = 1 — ft, i.e. the dynamics are the same 
for weak dissipation and short time, so it can be eas- 
ily understood that the geometric phases are nearly the 
same. When the dissipation is large, the items higher 
than the second-order can not be ignored for a quasi- 
periods. In this case, we set 8^0, i.e., the initial state 
vector in Bloch sphere is near the z axis. After a quasi- 
periods, the vector turns to the hemisphere containing 
the —z axis in Markovian case while it remains in the 
initial hemisphere in the non-Markovian one. This is the 
reason why the difference is so large between both pro- 
cesses with large dissipation rate at 8 — > 0, with which 



In this section, we will calculate the geometric phase 
of a two-level atom by using a memory kernel master 
equation with exponential memory phcnomcnologically. 
This equation may lead to a non-positive reduced den- 
sity matrix[23] and the positivity for a qubit has been 
discussed in Ref. [24j . In this paper, we focus on the ge- 
ometric phase only and do not take the positivity into 
account. 

We consider a two-level atom interacting with a vac- 
uum. An integrodifferential master equation containing 
the memory kernel in the interaction picture is[23l. [24| 

^-= I K(t-t')Cp(t')dt' = f K{t')£p(t-t')dt' (21) 
™ Jo Jo 

where C is the Liouvillian superoperator that may take 
the form 

£p = 7jJoC2p~ P& + - cr + (T~p - pa + (T~), (22) 

and K (f) is the memory kernel and here we choose an 
exponential memory phcnomcnologically 

K(t) = 7e~ 7t . (23) 

We call tr = — the memory time. 70 is the dissipation 
constant. One can solve this integrodifferential equation 
by taking its Laplace transform, determining the poles 
and inverting the solution in the standard way. The an- 
alytic solution of this equation [25( with initial condition 
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Eq.([2|) in Schrodinger representation is 



where £(R, r) is given by Eq.(|2il 
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sinh[-v/|l-4i?|] + cosh[-v/|l-4i?|] , 



(24) 



and i? = 22., r = 7*. This coefficients are valid for 4i? < 1 
and 2i? < 1. For AR > 1 and 2i? > 1, the form of these 
coefficients are acquired by substituting sinh[.] and cosh[.] 
with sin[.] and cos[.]. 

After the same procedure given above, we obtain the 
geometric phase acquired at time T = ^- as Eq. ([To"]) with 



2 °t 
cos — = 
2 



and 



(l + 7?-2£(E, t) cos 2 f) 2 



(1 + jj - 2£(i?, r) cos 2 f ) 2 + £ 2 (f , t) sin 2 9 



(25) 
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l-2cos 2 -£(i?,r 



e 2 (f ,r)sin 2 0. (26) 



The geometric phases of the system governed by the 
generalized Lindblad master equation as a function of 
9 are shown in Fig. ([3]), we plot this figure with expo- 
nential memory for different dissipation rate 70 and dif- 
ferent memory time 7. The geometric phase <&gp an d 
the polar angle 9 are plotted in units of %. We see 
from the figure that for weak dissipation(small dissipa- 
tion constant, 70 = 0.1), the three curves stand much 
close, even when the memory time is long(7 = 0.1). For 
large dissipation rate, long time memory lead to a large 
departure from the Markovian process. When the mem- 
ory time is shortened (7 = 10), the result tends to the 
Markovian. 



IV. GEOMETRIC PHASE WITH 
POST-MARKOVIAN MASTER EQUATION 

In Ref. [2y] , a post-Markovian master equation includ- 
ing bath memory effects via a phenomenological intro- 
duced memory kernel K (t) is derived. Compared with 
the memory kernel master equation, the advantage of 
post-Markovian master equation is that for a qubit it 
keeps the positivity of the density matrix for an expo- 
nential memory [25|. The general form of post-Markovian 
master equation for a two-level atom in zero-temperature 
reservoir is 




0.2 0.4 0.6 0.8 1 




(27) 



FIG. 3: Geometric phase of generalized Lindblad master 
equation with exponential memory as a function of 9. The 
geometric phase $gp and the polar angle 8 are plotted in 
units of 7r. The transition frequency u = 1 and time r = ^ZL. 
The parameter 70 is chosen: top: 70 = 0.1, bottom: 70 = 1. 
Solid line indicate 7 = 0.1, dash line 7=1, and dash-dot line 
7 = 10. 



where C is the Liouvillian superoperator and K(t) is the 
exponential memory kernel. It has been proved that 
the memory kernel master equation is a special case of 
post-Markovian master equation and one may derive the 
memory kernel master equation from the post-Markovian 
equation in the limit 70 <C 7[25|. 

The time dependent of the state for atom and the geo- 
metric phase after a quasi-periods have the same form as 



6 




"0 0.2 0.4 0.6 0.8 1 



9 

FIG. 4: The same as Fig. ([3} for post-Markovian master equa- 
tion. 



The geometric phase for the post-Markovian master 
equation is shown in Fig. ((4|. A similar feature as that 
from the memory kernel master equation can be seen 
in the figure. Furthermore when the dissipation rate is 
small, the figure is nearly the same as that obtained from 
the memory kernel master equation. The difference be- 
tween the two cases can be found in the large dissipa- 
tion case, e.g., the geometric phase departure from the 
Markovian process more largely with 7 = 0.1 than that 
in memory kernel master equation. 



V. CONCLUSION 

In summary, we have studied the geometric phases of 
a two-level system coupled to a non-Markovian environ- 
ment. The non-Markovian effects on the geometric phase 
are presented and discussed. We have chosen three differ- 
ent methods to describe the non-Markovian process. The 
common feature is that the geometric phase tends to zero 
with 9 approaches 71". As the memory kernel master equa- 
tion can be reproduced from the post-Markovian master 
equation in the limit of 70 <§C 7, the geometric phase cal- 
culated confirmed this point in this limit. Our results 
also show that for all the methods the non-Markovian 
effects change the geometric phase more strikingly when 
the dissipation rate and population of excited state in the 
initial state is large. This can be understood as the com- 
petition between the decoherence and the non-Markovian 
effects in the open system. 
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